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We analyze the behavior of systems of two and three qubits made by Josephson junctions, treated 
in the two level approximation, driven by a radiation mode in a cavity. The regime we consider is 
a strong coupling one recently experimentally reached for a single junction. Rabi oscillations are 
obtained with the frequency proportional to integer order Bessel functions in the limit of a large 
photon number, similarly to the case of the single qubit. A selection rule is derived for the appearance 
of Rabi oscillations. A quantum amplifier built with a large number of Josephson junctions in a 
cavity in the strong coupling regime is also described. 

PACS numbers: 74.50.+r, 03.65.Ud, 03.67.Lx, 42.50.Fx 
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Quantum computation 1 *^ promises a large improvement on the ability to execute demanding algorithms due to the 
large parallelism involved. Presently, it is not yet clearly understood how the hardware for a quantum computer should 
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ryj ' be realized. Several proposal relying e.g. on ion traps4 and NMR& have been devised but solid state devices seem to 
be very promising for this aim. Josephson junctions have been largely used both experimentally and theoretically for 
this goai§i 

In a couple of recent experiments Nakamura, Tsai and Pashkin? and Chiorescu et al£ were able to display the be- 
havior of a single Josephson junction strongly coupled to a radiation field and its suitability for quantum computation. 
The strong coupling regime has the advantage that can be theoretically described by a two-level modeli2*ii. 

Validity of the two-level model to describe the behavior of Josephson junctions has been proved by Stroud and Al- 
SaidiiSiiii. They considered the case of weak coupling where the so called rotating wave approximation does apply 14 . 
They also proved that, in some cases, a correction term to the model is needed while collapse and revival of Rabi 
■ oscillations is observed as also happens to similar models in quantum optics^i^. Similarities between solid state 
| devices and quantum optics systems are becoming increasingly meaningful^, 
—i ■ The model used by Stroud and Al-Saidi cannot be applied directly to the experiments of the Nakamura's group. 
Rather, we need to use the Dicke model— without any other approximations than the numbers of radiation modes 
and two-level systems. The first complete application of this kind has been given by us^S even if a study of this model 
in the strong coupling regime was started by Cohen- Tannoudji and his group 19 for a single radiation mode and a 
two-level atom. In this case all the terms in the interaction part of the Hamiltonian must be retained. 

In our paper— we were able to obtain the equations for the probability amplitudes in the strong coupling regime 
given the proper set of states. Two bands of levels were obtained and both intraband and intcrband transitions were 
shown. In the limit of a large number of photons in the cavity, the Rabi frequency is proportional to the integer 
number Bessel functions. We will recover this result below. Rather interestingly, the square of the amplitudes of 
I . the levels involved in the Rabi oscillations is a Poisson distribution. Rabi oscillations arise from the crossing of the 
energy levels of the Dicke model and appear between macroscopic superposition of charge states (known in the current 
literature as Schrodinger cat states). 

In the Dicke model, two-level systems are coupled by the radiation field. The same can happen with Josephson 
junctions 20 . So, our aim is to see how collective effects can emerge when more Josephson junctions are coupled by a 
cavity field. We will prove that Rabi oscillations emerge also in this case with the Rabi frequencies proportional to 
integer order Bessel functions. A selection rule arises for situations with more qubits as certain transitions are not 
allowed. In turn this means that the transition to some qubit configurations cannot be realized. Then, states involved 
in the Rabi oscillations are always macroscopic superpositions of entangled states between the radiation field, that 
can have a large number of photons, and the junctions. 

Finally, we extend the analysis to the thermodynamic limit showing that, in this case, a large number of Josephson 
junctions coupled by a cavity field in the strong coupling regime can be used to amplify vacuum fluctuations of 
the field to a large macroscopic field, assuming all the junctions in their ground state. Decoherence effects may be 
prominent in this case as we will discuss. This effect defines a quantum amplifier (QAMP) as proposed by us in the 
literaturo 14 i 21 i 22 . 

The paper is structured as follows. In sec[n]we introduce the model describing a number of Josephson junctions 
interacting through a cavity field. In sec IIIII we analyze explicitly the cases for one, two and three coupled qubits 
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obtaining the Rabi frequencies and proving, in the limit of a large number of photons, the proportionality with integer 
order Bessel functions. In sec II VI we discuss the thermodynamic limit, that is we see the physics of a large number of 
Josephson junctions in a cavity proving that an amplification of the vacuum fluctuations of the radiation field in the 
cavity is obtained. This new device we term QAMP. Finally, in sec0 conclusions are given. 



II. DESCRIPTION OF THE MODEL 



The model we should consider for a single Josephson junction, treated in the two-level approximation, in a cavity 
field is (here and in the following % = 1) 

H = AS Z + ua)a + 2gS x {a ] + a) (1) 

being A the separation between the ground and the first excited state in a Cooper paisi2*A£, ui the frequency of the 
field in the cavity, g the coupling between the junction and the cavity, a' and a creation and annihilation operators, 
S x , S z = % con a x , a z Pauli matrices. 

Eq.Q is just the Dicke model that we have specialized to a single qubit with Pauli matrices. This generalizes 
immediately to any number of qubits putting 

1 N 

t=l 

1 N 

S z = -j^Z^z.i 

i=l 

being N the number of qubits. 

One may think that eq.Q differs from the Hamiltonian currently used to describe Josephson junction qubits in 
experiments^ but these Hamiltonians are equivalent as we can interchange a x and a z with the unitary transformation 
e l ~°y and what one gets, in the analysis in the strong coupling regime, is a shift to the energy levels. So, both 
Hamiltonians give rise to the same physics and we can safely work with the Dicke model as currently known. 

In order to work out the analysis of the model in the strong coupling regime proper to the experiments carried out 
so far, we follow the same approach used in RefiiS. We split the Hamiltonian JIJ as 

Hq = ua'a + 2gS x (<r + a) (3) 
Ht = AS Z 

being Hq the unperturbed Hamiltonian and H\ the perturbation. The unperturbed part Hq can be immediately 
diagonalized by the eigenstates, as can be directly verified treating S x as a c-number, 

\[n};S,S x )^e^ a ' a ^\n)\S,S x ) (4) 

being S x the component of the spin along x axis chosen as the axis of quantization, \S, S x ) the corresponding eigenstates 
and a^a\n) = n\n) and n = 0, 1, 2, . . . an integer. Energy eigenvalues are given by 

E„. Sx =nu (5) 

u> 

and we observe a degeneracy between positive and negative values of the spin components. Then, we can introduce 
the unitary evolution operator 

( f«i£iV 

^ — , — ? nu — t 

M*)=5> V ' \[n];S,S x ){[n};S,S x \. (6) 

We can apply this unitary operator to the initial Hamiltonian Q and the problem is reduced to the Hamiltonian 

— _ _ — _ i \ muj — t —% I nuj — t 

#=]T5>V J e V - J A\[m];S,S x }([n];S,S' x \([m];S,S x \S z \[n];S,S' x ). (7) 

m,S x n.S' 
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The matrix elements can be evaluated by noting that S z — + 2 - and one gets 



([m];S,S x \S z \[n];S,S' x ) = i (8s a ,s>.+iVS(S + 1) - S> X (S> X + l)Af-„ (8) 
- ^Si-iv 7 ^ + 1) - ^ - l)M+„) 

where we have puliiS 

ML = M^<"V> = Me^ (±*°Y™ (9) 

V to! V uj I \ u> z I 



being L™ " Laguerre polynomials. 

Now the Hamiltonian (JJJ can be split in two parts to give a diagonal part 

H = ]T ^AM„ (e-'^'^WSiS + 1) - S X (S X - l)\[n];S, S x )([n];S, S x - 1| (10) 



2i 

,. 4g 2 (2S :c +l) 

e 



VS , (5+l)-5 a: (5 a : + l)|[n];5, S x > (N ; S, 5, + 1 1) , 



having set 



and an off-diagonal part 



M„ = M*„ = e-^L„ ( -Jj. ) , (11) 



ffi = E ^ e " i( "" m)ut (e-' 4a2<2 ' x ' 1) t M m „ v /^(g+ 1) - S X (S X - l)\[m];S,S x )([n];S,S x - 1| (12) 



■ 4 fl 2 (2S a; + l) 



e 4 ^'M+„ V^(5 + 1) - 5,(5, + 1)| [to]; 5, &><[«]; S, 5 a + 1| . 



In order to obtain the equations for the probability amplitudes, we diagonalize Hamiltonian (|10f> obtaining the 
energy eigenvalues, the eigenstates and the geometrical part of the phase as we can have time-dependent eigenstates, 
in agreement to the general approach outlined in Rcf. 10 23 . Then, by interaction picture, the amplitude equations are 
obtained using Hamiltonian Q12[l. 

In the analysis we pursue in the following we will see that Hamiltonian l|10|) has eigenvalues s x AM n having s x the 
same values of the spin projection S x and one can have transitions between the eigenstates only if the selection rule 
As, = 0, ±1 holds. We recognize that each level s x develops a band of levels numbered by integer n, so the selection 
rule provides intraband (s x — s x ) and interband (s x ^ s x ) permitted transitions at resonance respectively. 

In the limit of a large number of photons we can show, in all cases below, that the Rabi frequencies are proportional 
to integer order Bessel functions. This can be accomplished using the relation^ 

J a (2V^) =e' f (^)~ L a n {x) (13) 
that holds in the limit of n going to infinity. 



III. ANALYSIS OF THE MODEL 



We apply the procedure outlined in the previous section to one, two and three qubits to see in details the physics 
of the model. The one qubit case has been already discussed in literature^ but we present it here again to show the 
way our method works to make this paper self-contained. 
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A. One qubit 

In the one qubit case we have S x = being S = |. The diagonal part of the Hamiltonian wil be 



r 1 1 1 



r l 1 1 



11 

• n]; 2'2 



(14) 



that is time-independent in this case. Similarly, we have the off-diagonal part 



1 1 

2' 2 



1 1 



- M+ 

mn 



l l 



m;- 



ra ; 



1 1 

2' 2 



(15) 



This operator is Hcrmitian as can be verified using the fact that sum indexes are dummy. 

Hamiltonian l|14|l can be immediately diagonalized with eigenvalues -E n ,s x = s x AM n being s x = |, — i corresponding 
to the eigenstates 



\n;s x ) 



V2 ( 



1 1 



i 



1 1 

2' 2 



(16) 



that are not eigenstates of the spin projection S x . Then, we can write the equation for the probability amplitudes by 
looking for a time-dependent solution in the form 



i^)) = E e " 



■iE k 



t Ck,s x (t)\k;s x ) 



(17) 



where no contribution enters due to geometrical phases as Hamiltonian 114fl is time-independent. The Schrodinger 
equation takes the form 



A 



-i\E„ 



(18) 



We recognize the resonance conditions 



E m ,s x + (n - m)uj = 



(19) 



with As x = for intraband transitions and As x = ±1 for interband transitions as said above. These resonance 
conditions correspond to crossing of the energy levels of the Dicke model. Energy levels are given bj*^4 



E,, 



is 2 a 2 

UUJ — + St A M„ 



(20) 



where we can recognize the band structure and the degeneracy for s x = ±i that is removed by the last term. Then, 
resonance conditions can be straightforwardly interpreted as crossings of energy levels. Such crossings produce Rabi 
oscillations as experimentally observed. 

Finally, let us check the Rabi frequencies for the single qubit case. It is not difficult to realize that we have for 
intraband transitions As x = from ca. (|18|) . assuming the rotating wave approximation due to the resonant terms, 



(M+„ + M- n ) = A (m 



and for interband transitions As T = ±1 



H' x = -A (M+ n - M~ 



cosh 



sinh 



2.9 



(a -a^ 



— (a — a 1 ) 



(21) 



(22) 



Now, we can use eq. (|!|J to show that transitions with states differing by an odd number of photons, m = n + 2N + 1 are 
permitted for intraband transitions and even number of photons are involved m = n+ 2N in interband transitions^. 
Using Sterling formula and ea. 1)13(1 one gets, in the limit of Fock states with a large number of photons, 



IZi w AJ 2 n+i 



w ) 



(23) 
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and 

K\ « AJ 2Ar (^) . (24) 

completely in agreement with experimental results by Nakamura, Tsai and PashkirA As also experimentally observed, 
crossing at m = n can happen belonging to interband transitions with TV = 0. Further details for a single qubit case 
are given in Ref^S. At present, this is the only case where a comparison with experiment can be realized. 



B. Two qubits 

For the case of two qubits one has 



H = ^i^AM„(e- <4 #*|[n];l,l)([n];l,0| (25) 



e 



[n];l,0)([n];l,-l| 



- e <1 S- t |[n];l J 0><[n];l J l| 

- e-^IHl.-lXMjl.Ol 

Now, the Hamiltonian shows time-dependence and it is not straightforward to interpret the resonance conditions as 
energy level crossings. On the same ground one has 

Hi = e-«- m ^V2A (e-^M-Jim]; 1, 1) ([n]; 1,0] (26) 

+ e^M^Km]; 1,0) <[n];l,-l| 

- M+ n \[m}; 1,0) ([n]-l,l\ 

- e^M+JH;!,-!)^];!^! 

One can immediately diagonalize Hamiltonian (|25|l obtaining eigenvalues E njSx — s x AM n with s x = 1,0,-1. The 
corresponding eigenvectors are 

\n;l,t) = l e -^*|[n];l,l) + i=|[n];l,0)-i e -^*|[n];l,-l) (27) 

|n;0,t) = ^|M;l,l) + i=|[n];l,-l) 

\n;-l,t) = -I e -^|[n];l,l) + i=|[n]; 1,0) + i e -^*|[n]; 1,-1) 

that do not give any geometrical phase as (n;l,t\idt\n;l,t) = (n;0,t\idt\n;0,t) = (n; — 1, i|i$t|n; — 1, t) = 0. So, we 
can put again 

W*)> = I>-^-'- t <*,.. w> (28) 

and obtain the equations for the amplitudes 

A 



ic k 



1 (t) = Y /T (M^ n +M+Je-^' E ^ + ^ k ^ t c n At)+ (29) 



ic k , (t) = J2^Ti {M ^ ~ M L>- ilEk ^ +(n - k)U]t Cn^ 



i[E h ,- 1 + =L-+(n-k)u]t 



n,-l(*) 



C 

iCk.-iW = £-j( M *n + M+)e-^— £ — + ("^] t c„,_ 1 (t) + 

n 

^2^(M- n - M+ )e < ^.- I +^-("- fe H t c Tl ,o(t) 

where we can recognize again the effect of the selection rules As^ = 0, ±1 on the permitted transitions, originating 
from the resonance conditions 

E n ,±i-Ek,±i + (n-k)u = (30) 
4<? 2 , 

Ek,±i H (n — s)td = 

for intraband and interband transitions respectively, giving rise to Rabi oscillations in the two-qubit system. Again, 
we can interpret this resonance conditions as originating from the crossings of the energy levels for the Dicke model 

As 2 a 2 

E n ^ =noj- —22- + s x AM n (31) 

CO 

that is degenerate with respect to s x = ±1 but the degeneracy is removed by the last term. The Rabi frequencies 
can be computed from eas. 1)29(1 with the rotating wave approximation imposing the resonance conditions 1)30(1 and are 
given by, for intraband transitions 

n 2 = n t (32) 

and for interband transitions 

n' 2 = VZR/i. (33) 

These frequencies, in the limit of a large number of photons, displays dependence on the integer order Bessel functions, 
being proportional to IZi and IZ'i like in the single qubit case. 

A cooperative effect arises from the coherent behavior of both the junctions entering entangled states and producing 
Rabi oscillations. 

So, the main conclusion of this section is that the two qubits act collectively producing Rabi oscillations, exactly 
as in the case of a single qubit. Such oscillations appear between entangled states to be considered macroscopic. 



C. Three qubits 

In order to gain further insight into the physics of Josephson junctions in this case, we analyze a system with 
three qubits. This case is rather different from the preceding ones as Hq depends on time and we have to apply 
a theorem for strong coupling that imposes formally the adiabatic theorem at the leading order™. So there is a 
geometrical contribution to the phases. This makes no straightforward to interpret the resonance conditions as energy 
level crossings. We have 



Ho = ^ AM « ( e- i£ ^V3 



2i 



3 3 
W; 2'2 



3 1 

n]; 2'2 



(34) 



3 1\ /. . 3 1 
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and 



+ 2M~ 



3 3 

2' 2 



[m] ; I , - > / fn 1 



'V3M: 



2' 2/ \ L '" J ' 2' 2 

3 1\ /, , 3 3 
n 



2' 2 
3 I 

2'^2 



3 1 



^ ; 2'-2/\ LnJ; 2'2 



3 1 



e- i *-ir t V3M+ 



3 3 



[mJ; 2'"2/\ [71J; 2'"2 



2' 2 
. 3 3 
l]; 2'~2 



3 1 



It is straightforward to diagonalize Hamiltonian (|34|l with the eigenstates 



3 1 

2' 2 



(35) 



|n;s x ,t) = j3(s x ) 



+ i \s x 



-ze 



4sx 



2s T 



3 3 
2' 2 



3 1 

2' 2 



3 1\ . l5 £! t V3, , 
^2'-2/- 6 " "I 1 4,? 



3 3 

1 J ' 2' 2 



(36) 



being 



/3(sx) 



(37) 



and s x = |, h, — h, — §. As in the other cases, we get entangled macroscopic states. In this case we have geometrical 
phases as given by 



2p 



2 3 



3 , 27 
2sT ' 16sl 



W + _3 



a: ^ 16sJ ^ 64s 



27 , 1 

i I A 



(38) 



that reduces to — for s x — §, — 1 and — for s x = i, —I. 
For the given Hq and Hi we seek a solution like 



(39) 



yielding the equations for the amplitudes 



,(*) = £■ 



,-t[£n,, a> -^,S a ,+(n-m)u>+7a a -7« a ,]< 



A0(s x )/3(s x ) [a(s x ,s a )M+ n + a(s :c ,s x )M- n ] c„, s ,(i) (40) 



having set 



3_ 3 / J3_\ / J_ 

4 1 4^ Sl 4,s K 



(41) 



It is not difficult to verify, also in this case, that the selection rules As x = 0,±1 hold corresponding to intraband 
(s x = s x ) and interband (s x 7^ s x ) transitions respectively with the resonance conditions 



E n ,s x - E m> s x +(n— m)uj + j Sx - j Sal = 0. 



(42) 
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for intraband transitions and 



The Rabi frequencies can be computed from eas. (|40() with the rotating wave approximation imposing the resonance 
conditions l|42l) and can be given explicitly as follows 

3 

72-3/2,3/2 = 72.-3/2,-3/2 = ^1 (43) 

1 

72-1/2,1/2 — 72—1/2,-1/2 — 2^1 

"^•3/2,1/2 = 72 1/2 ,3/2 = 

72-3/2,-1/2 = 72_i/ 2 ,-3/2 = -^-T^-'l 
72-1/2,-1/2 = 72._l/2,l/2 = 72. 1 

for interband transitions. Also for three qubit the Rabi frequencies are proportional to 72 1 and 72.' i , so in the limit 
of a large number of photons one has that these frequencies are given by integer order Bessel functions. This rule 
applies to all cases. 

As also seen for one and two qubits we can see again a collective effect. The behavior of the three qubits is coherent 
producing entangled states with Rabi oscillations. Apart from numerical factors, the Rabi frequencies are the same 
as in the single qubit case. It would be interesting to extend the analysis to higher spins to sec if this is a rule. 

IV. THERMODYNAMIC LIMIT: QUANTUM AMPLIFIER 

To analyze this case, we limit the study to the leading order as higher orders become increasingly less important 
as the number of Josephson junctions increases^. This permits us to write down immediately the unitary evolution 
being given by Hq in eq.®- It is important to emphasize that, in order to unveil this effect, we should avoid random 
phases putting all the junctions in the same state. This makes the case somewhat different from above but shows 
again a collective effect possible to be observed. 

Now, let us assume that the field in the cavity is in the ground state and all the Josephson junctions are in their 
lower state externally imposed. Unitary evolution gives 



( g 2 iV 2 V 

i 1 nuj— — I i 



N N\ »V [Nq\ n 1 , 



n 

So, leaving aside the contributions of the Josephson junctions, one is left with 

Ngr t\ x — a —i(nw— 9 " )t N 2 g 2 ( Na\ n 1 

that is, a coherent state with a parameter proportional to N. In the large N limit this represents a classical state of 
the radiation field 2 ^ and we see that we have amplified quantum fluctuations of the ground state of the radiation field 
to a classical level. The interesting aspect of this already known result is that it can be practically realized through 
Josephson junctions that lend themselves to realize this kind of classical state out of a fully quantum initial state. 
One realizes a QAMP (Quantum Amplifier) 14 ! 22 . 

The device can be realized without too much care about taking just the limit N — > oo or N — > oo, V — ► oo, being 
V the volume of the cavity, and N/V — constant as we are in the strong coupling limit. 

This represents a very peculiar collective effect of a large number of Josephson junctions in a cavity in the strong 
coupling regime, the same regime devised in the experiments of Nakamura's group 6 '- for a single qubit. 

V. DISCUSSION AND CONCLUSIONS 



The above discussed collective effects can have their limit in the appearance of decoherence. The experiments of 
Nakamura's group display the decay of Rabi oscillation but the decay time can be enough long to permit quantum 
computation. Neither the nature of this decoherence effects nor the way it scales with the number of Josephson 
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junctions are known as far as we know. It should be expected that the same should happen for the QAMP producing 
a classical radiation field. The origin of decohcrcncc on Joscphson junctions is to be understood but it could not 
be excluded that revival and collapse effects of Rabi oscillations may be at work as happens in the weak coupling 
regim o i . 

We have shown how, in the strong coupling limit as devised by Nakamura's group, several collective effects as 
Rabi oscillations and quantum amplification could appear in Josephson junctions coupled by a radiation field. Rabi 
oscillations, in the limit of a large number of photons, have frequencies proportional to integer order Bessel functions 
and a selection rule applies limiting transitions. This effects could be useful for applications in quantum computation 
or to use a large number of Josephson junctions to produce a laser field. 



* marcofrasca@mclink.it 
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